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Abstract. 

Motivated by fundamental questions about the loss of phase coherence 
at low temperature we consider relaxation, dephasing and renormalization 
effects in quantum two-level systems which are coupled to a dissipative 
environment. We observe that experimental conditions, e.g., details of the 
initial state preparation, determine to which extent the environment leads 
to dephasing or to renormalization effects. We analyze an exactly solvable 
limit where the relation between both can be demonstrated explicitly. We 
also study the effects of dephasing and renormalization on response func- 
tions. 



1. Introduction 

The dynamics of quantum two-level systems has always been at the focus of 
interest, but recently attracted increased attention because of the prospects 
of quantum state engineering and related low-temperature experiments. A 
crucial requirement for many of these concepts is the preservation of phase 
coherence in the presence of a noisy environment. Typically one lacks a 
detailed microscopic description of the noise source, but frequently it is 
sufficient to model the environment by a bath of harmonic oscillators, with 
frequency spectrum adjusted to reproduce the observed power spectrum. 
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The resulting 'spin-boson' models have been studied in the literature, in 
particular the one with bilinear coupling and 'Ohmic' spectrum, but spin- 
boson models with different power spectra appear equally interesting in 
view of several experiments. In this article we, therefore, analyze spin-boson 
model with general power spectra with respect to relaxation, dephasing and 
renormalization processes at low temperatures. We study the dephasing of 
nonequilibrium initial states. We show how the state preparation affects 
the effective high-frequency cut-off and separates renormalization from de- 
phasing effects. We also demonstrate how dephasing processes influence 
low-temperature linear response and correlation functions. 

2. Spin-boson model 

In this section we review the theory and properties of spin-boson models, 
which have been studied extensively before (see the reviews [4, 11]). A 
quantum two-level system is modeled by a spin degree of freedom in a 
magnetic field. It is coupled linearly to an oscillator bath representing the 
environment. The total Hamiltonian is 

U = Wctrl + Oz E C i(°j + 4) + Hh ' W 

3 

where the controlled part is 7Y c tri = ~\B Z (J z~\B x a x = — ^AE (cos 6 a z + 
sin 6 a x ), the oscillator bath is described by Ti^, = J2j a } a j > an d the 

bath 'force' operator X = J2j c j( a j + a }) ^ s assumed to couple linearly to 
a z . For later convenience H c t r i has also been written in terms of a mixing 
angle 9 = tan~ 1 (B x / 'B z ), depending on the direction of the magnetic field, 
and the energy splitting of the eigenstates, AE = yj B'%. + . 

In thermal equilibrium the Fourier transform of the symmetrized corre- 
lation function of the force operator is given by 

S x {u) ee ([X(t),Xtf)] + ) u = 2ftJMcoth^ . (2) 

Here the bath spectral density has been introduced, defined by J(co) = 
f J2j c ] — ujj) . At low frequencies it typically has a power-law form up 
to a high-frequency cut-off uj c , 

J(u) = \ha w£- a w a e(w c - uj) . (3) 

Generally one distinguishes Ohmic (s = 1), sub-Ohmic (s < 1), and super- 
Ohmic (s > 1) spectra. In Eq. (3) an additional frequency scale ujq has been 
introduced. Since it appears only in the combination acv^ s it is arbitrary 
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and in Ref. [4] has been chosen equal to the high-frequency cut-off. Here 
we prefer to distinguish both frequencies, since the cut-off ui c will play an 
important role in what follows. 

The spin-boson model has been studied mostly for the case where an 
Ohmic bath is coupled linearly to the spin. One reason is that linear damp- 
ing, proportional to the velocity, is encountered frequently in real systems. 
Another is that suitable systems with Ohmic damping show a quantum 
phase transition at a critical strength of the dissipation, a CT ~ 1. On the 
other hand, in the context of quantum-state engineering we should concen- 
trate on systems with weak damping, but allow for general spectra of the 
fluctuations. 

Spin-boson models with sub-Ohmic damping (0 < s < 1) have been 
considered earlier [4, 11] but did not attract much attention. It was argued 
that sub-Ohmic dissipation would totally suppress coherence, transitions 
between the states of the two-level system would happen only at finite 
temperatures and would be incoherent. At zero temperature the system 
should be localized in one of the eigenstates of a z , since the bath renor- 
malizes the off-diagonal part of the Hamiltonian B x to zero. This scenario 
is correct for intermediate to strong damping. It is, however, not correct 
for weak damping. Indeed the 'NIBA' approximation developed in Ref. [4] 
fails in the weak-coupling limit for transverse noise, while a more accurate 
renormalization procedure [3] predicts damped coherent behavior. In the 
context of quantum-state engineering we are interested in precisely this 
coherent sub-Ohmic regime. We will demonstrate a simple criterion which 
allows to define a border between coherent and incoherent regimes. 

A further reason to study the sub-Ohmic case is that it allows us to 
mimic the universally observed 1/ / noise. For instance, a bath with s = 
and J (a;) = (ir/2)ah~oJo produces at low frequencies, Tiu <C &bT, a 1// noise 
spectrum Sx{oj) = Ef^/\uj\ with Ef^ = l-Kcthuok^T . Since frequently 
nonequilibrium sources are responsible for the 1// noise, the temperature 
here should be regarded as a fitting parameter rather than a thermodynamic 
quantity. 

Below we will also consider the super-Ohmic case (s > 1) as it allows 
us to study renormalization effects clearly. 

3. Relaxation and dephasing 

We first consider the Ohmic case in the weak damping regime, a<l. Still 
we distinguish two regimes: the 'Hamiltonian-dominated' regime, which is 
realized when AE is large enough, and the 'noise-dominated' regime, which 
is realized, e.g., at degeneracy points where AE — > 0. The exact border 
between both regimes will be specified below. 
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In the Hamiltonian-dominated regime it is natural to describe the evo- 
lution of the system in the eigenstates of H c tr\, which are 

|0)= cos °- |T) + sin 6 -\[) and |1) = - sin °- \\) + cos °- ||) . (4) 

Denoting by t x and t z the Pauli matrices in the eigenbasis, we have 

H = -^AEt z + (sin0 t x + cos6> r z ) X + H h . (5) 



Two different time scales describe the evolution in the spin-boson model. 
The first, the dephasing time scale t v , characterizes the decay of the off- 
diagonal elements of the qubit's reduced density matrix p(t) in the preferred 
eigenbasis (4), or, equivalently of the expectation values of the operators 
t± = (\/2)(T x ±iT y ). Frequently dephasing processes lead to an exponential 
long-time dependence, 

<T±(t)> = tr [r±p(t)] oc (r±(0)) e* iAEt / h e"*^ , (6) 

but other decay laws occur as well and will be discussed below. The second, 
the relaxation time scale r re i ax , characterizes how diagonal entries tend to 
their equilibrium values, 

<r z (t)> - <r 2 (oo)> oce-*/^'- , (7) 

where (r 2 (oo)) = tanh(A£/2/c B T). 

In Refs. [4, 11] the dephasing and relaxation times were evaluated in a 
path-integral technique. In the regime a <C 1 it is easier to employ the per- 
turbative (diagrammatic) technique developed in Ref. [9] and the standard 
Bloch-Redfield approximation. The rates are 

r re i ax = r-^ = ^ sin 2 6S X {u = AE/h) , (8) 

= t- 1 = \ r rclax + cos 2 es x (u> = o) , (9) 

where Sx{w) = 7rah 2 uj coth(hu>/2kBT). We observe that only the 'trans- 
verse' part of the fluctuating field X, coupling to t x and proportional to 
sin 6, induces transitions between the eigenstates (4) of the unperturbed 
system. Thus the relaxation rate 1 (8) is proportional to sin 2 9. The 'longi- 
tudinal' part of the coupling of X to t z , which is proportional to cos 9, does 

x ThG equilibration is due to two processes, excitation |0) — ► |1) and relaxation |1) — ► 
|0), with rates oc {X(t)X(t')} w= ±A. E /h- Both rates are related by a detailed balance 

condition, and the equilibrium value (r z (oo)) depends on both. On the other hand, T^iax 
is determined by the sum of both rates, i.e., the symmetrized noise power spectrum Sx ■ 
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not induce relaxation processes. It does, however, contribute to dephasing 
since it leads to fluctuations of the eigenenergies and, thus, to a random rel- 
ative phase between the two eigenstates. This is the origin of the 'pure' de- 
phasing contribution to Eq. (9), which is proportional to cos 2 9. We rewrite 
Eq. (9) as T v = ±r re i ax + cos 2 6T* , where T* = S x (w = 0)/fi 2 = 2irak B T/h 
is the pure dephasing rate. 

The pure dephasing rate T* characterizes the strength of the dissipative 
part of the Hamiltonian, while AE characterizes the coherent part. The 
Hamiltonian-dominated regime is realized when AE » r* while the noise- 
dominated regime is realized in the opposite case. 

In the noise-dominated regime, AE <C T* , the coupling to the bath is 
the dominant part of the total Hamiltonian. Therefore, it is more convenient 
to discuss the problem in the eigenbasis of the observable a z to which 
the bath is coupled. The spin can tunnel incoherently between the two 
eigenstates of a z . One can again employ the perturbative analysis [9], but 
use directly the Markov instead of the Bloch-Redfield approximation. The 
resulting rates are given by 

Treiax = Bl/Y% = B 2 J(2irhak B T) (10) 
r v = T* = 2-irhak B T . 

In this regime the dephasing is much faster than the relaxation. In fact, as 
a function of the coupling strength a the dephasing and relaxation rates 
evolve in opposite directions. The a-dependence of the relaxation rate is 
an indication of the Zeno (watchdog) effect [2] : the environment frequently 
'observes' the state of the spin, thus preventing it from tunneling. 

4. Longitudinal coupling, exact solution for factorized initial con- 
ditions 

The last forms of Eqs. (8) and (9) express the two rates in terms of the noise 
power spectrum at the relevant frequencies. These are the level spacing 
of the two-state system and zero frequency, respectively. The expressions 
apply in the weak-coupling limit for spectra which are regular at these 
frequencies. For the relaxation rate (8) the generalization to sub- and super- 
Ohmic cases merely requires substituting the relevant Sx{w = AE/U). 
However, this does not work for the pure dephasing contribution. Indeed 
Sx{w = 0) is infinite in the sub-Ohmic regime, while it vanishes for the 
super-Ohmic case and the Ohmic case at T = 0. As we will see this does 
not imply infinitely fast or slow dephasing in these cases. To analyze these 
cases we study the exactly soluble model of longitudinal coupling, 9 = 0. 

Dephasing processes are contained in the time evolution of the quantity 
(r+(i)) obtained after tracing out the bath. This quantity can be evaluated 
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analytically for 9 = (when r + = a + ) for an initial state which is described 
by a factorized density matrix p(t = 0) = /5 S pin ® Pbath (which implies 
that the two- level system and bath are initially disentangled). In this case 
the Hamiltonian H = — ^AEa z + a z X + Tib is diagonalized by a unitary 
transformation by U = exp (— icr z $/2) where the bath operator $ is defined 
as 

• = «Es£<<4-«,>- (») 



The 'polaron transformation' yields 7i = UTCU 



lAEa z + H h . It 



has a clear physical meaning: the oscillators are shifted in a direction de- 
pending on the state of the spin. Next we observe that the operator cr + is 
transformed as 5+ = XJo+U^ 1 = e _l *<7+ , and the observable of interest 
can be expressed as 



<a+(t)) = Tr [p(t = 0)o+(t)] = Tr \Up(t = 0)U- 1 e~ i ^a + 



(12) 



The time evolution of <&(£) = e* Wb *$e _iWb * is governed by the bare bath 
Hamiltonian. After some algebra, using the fact that the initial density 
matrix is factorized, we obtain (o+(t)) = V(t) e~ iAEt <ct+(0)) where 



7>(t) = Tr e 



a i*(0)/2. 



(13) 



The expression (13) applies for any initial state of the bath as long as 
it is factorized from the spin. In particular, we can assume that the spin 
was initially (for t < 0) kept in the state ||) and the bath had relaxed 
to the thermal equilibrium distribution for this spin value: /Sbath = Pi = 



Z^e-^T, where U 



Wb + EjCj(o- + a}) 

- l ®/ 2 } where the density matrix of the 
decoupled bath is given by pb = Z h ~ 1 e~ f3 ' Hh , and the function V(t) reduces 
to 

Vit) P(t) = Tr (e-**^ p b ) . (14) 

This expression (with Fourier transform P{E)) has been studied extensively 
in the literature [4, 8, 7, 6, 1]. It can be expressed as P(t) = exp K (t), where 



In this case we can rewrite the 



density matrix as Pbath = e l */ 2 /?t>e 



Kit) 



coth 



\2k B T 



(cos cot — 1) — i sin tot 



(15) 



For an Ohmic bath at non-zero temperature and not too short times, 
t > h/k B T, it reduces to ReK(t) w -S x (uj = 0)t/h 2 = -2irak B Tt/h, 
consistent with Eq. (9) in the limit 9 = 0. 
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For 1/uj c < t < h/k^T, and thus for all times at T = 0, one still finds 
a decay of (a+(t)) governed by Re K(t) «-2a ln(u; c i), which implies a 
power-law decay 

(a + (t)) = (co c t)-^e-* AEt / h (v + (0)) ■ (16) 

Thus even at T = 0, when Sx{w = 0) = 0, the off-diagonal elements of the 
density matrix decay in time. All oscillators up to the high-frequency cut- 
off uj c contribute to this decay. The physical meaning of this result will be 
discussed later. We can also define a cross-over temperature T* below which 
the power-law decay dominates over the subsequent exponential one. A cri- 
terion is that at t = fijk^T* the short-time power-law decay has reduced 
the off-diagonal components already substantially, i.e {<j+{Ti/k-B,T*)) = 1/e. 
This happens at the temperature kT* = fru) c exp(—l/2a). Thus the dephas- 
ing rate is T* = k^T* jfi for T < T* and T* = lixak^T jU for aT > T* . 

For a sub-Ohmic bath with < s < 1 due to the high density of 
low-frequency oscillators exponential dephasing is observed even for short 
times, \(a + (t))\ oc exp[— a(woi) 1_s ] for t < fi/k^T, as well as for longer 
times, |(<7+(i))| oc exp[— a Tt (wot) 1 " 5 ] for t > hjk-^T. In the exponents of 
the short- and long-time decay laws we have omitted factors which are of 
order one, except if s is close to either 1 or 0, in which case a more careful 
treatment is required. The dephasing rates resulting from the decay laws are 
r* oc T* = a 1 /( 1 -*) Wo (cf. Ref. [10]) for T <T* and r* oc (aT/w ) 1/{2 ~ s) ^o 
for T > T*. Again the crossover temperature T* marks the boundary be- 
tween the regimes where either the initial, temperature-independent decay 
or the subsequent decay at t > h/k-^T is more important. 

These results allow us to further clarify the question of coherent vs. 
incoherent behavior in the sub-Ohmic regime. It is known from earlier 
work [4, 11] that the dynamics of the spin is incoherent even at T = 
if AE <C o}-/^-~ s ^7ujJq (only the transverse case 6 = ir/2 was considered). 
This condition implies AE <C r* i.e., it marks the noise-dominated regime. 
We are mostly interested in the opposite, Hamiltonian-dominated regime, 
AE 3> hT*. The NIB A approximation used in Ref. [4] fails in this limit, 
while a more accurate RG study [3] predicts damped coherent behavior. 

Finally we discuss the super-Ohmic regime, s > 1. In this case, within 
a short time of order w" 1 , the exponent Re Kit) increases to a finite value 
Re K(co) = —a{uj c /uJo) s ~ l and then remains constant for a;" 1 < t < fi/k^T 
(we again omit factors of order one and note that the limit s — > 1 requires 
more care). This implies an initial reduction of the off-diagonal element 
|(<7+(i))| followed by a saturation at \(a + (t)}\ oc exp[— a(a; c /wo) s_1 ]- For 
t > fi/k^T an exponential decay develops, |(<r+(£))| oc exp[— aTt (woi) 1_s ], 
but only if s < 2. This decay is always dominant and, thus, there is no 
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TABLE 1. Decay of |(<7+(t))| in time for different bath spectra. 
Ohmic: J(lo) = ^aujQ(ujc — to): 

\<o + (t)>\ 



for to' 1 < t < T- 1 |(a+(t)}| » (co c t)- 2a 
for t>T- x |(a+(t))|^e- 2 ' rQTt 



power 
law 



exp 



Sub-Ohmic: J(u) = § auj s w s 6(we - w) and < s < 1: 



for W" 1 < t < T" 1 |(ff+(t))| ^e- a( " ot)1 3 
for ^T- 1 |(<T+(t)>| » e-^tf-ot) 1 - 



1/T 



l«y + (t)>i 



exp 



exp 



Super-Ohmic: J(w) = fawj s cj s 9(lj c - w) and 1 < s < 2: 



l/T 



l«y + (t)>| 



forw- 1 <t<T- 1 |(a+(t))| » e -"(^/-o) s 1 
for ^T- 1 Ktr+WJI^e- 05,1 ^') 1 " 



partial 



• Soft gap: J(w) = fa^o ^'©(tdc — <-o) and s > 2: 



l/T 



l«y + (t)>i 



for to' 1 < t 



KM«)>I 



-ci(u c /^ ) s 1 



partial 



l/T 



crossover in this case, i.e., T* = 0. For s > 2 there is almost no additional 
decay. 

The results obtained for different bath spectra are summarized in Ta- 
ble 1. 



DEPHASING AND RENORMALIZATION 



9 



5. Preparation effects 

In the previous section we considered specific initial conditions with a fac- 
torized density matrix. The bath was prepared in an equilibrium state char- 
acterized by temperature T, while the spin state was arbitrary. The state 
of the total system is, thus, a nonequilibrium one and dephasing is to be 
expected even at zero bath temperature. In this section we investigate what 
initial conditions may arise in real experiments. 

The fully factorized initial state just described can in principle be pre- 
pared by the following Gedanken experiment: The spin is forced, e.g., by a 
strong external field, to be in a fixed state, say ||). The bath, which is kept 
coupled to the spin, relaxes to the equilibrium state of the Hamiltonian H^, 
e.g., at T = to the ground state |g|) of H\. Then, at t = 0, a sudden pulse 
of the external field is applied to change the spin state, e.g., to a superposi- 
tion -^(||) + |J,)). Since the bath has no time to respond, the resulting state 

is |i) = ^ (|t) + |D) (g) |g|). Both components of this initial wave function 
now evolve in time according to the Hamiltonian (1). The first, ||) (g) |g|), 
which is an eigenstate of (1), acquires only a trivial phase factor. The time 
evolution of the second component is more involved. Up to a phase factor 
it is given by ||) (g> exp(-iHit/h) |gj) where Hi = Tib - J2j c ji a j + a ])- 
As the state |g-f) is not an eigenstate of entanglement between the 
spin and the bath develops, and the coherence between the components 
of the spin's state is reduced by the factor | (g|| exp(— iTCit/h) |gj) | = 
| (g |exp(-^(t))exp(i$)|g ) | = \P Uc (t,T = 0)| < 1. The function P(t) 
was defined in Eq. (14). The subscript to c is added to indicate the value of 
the high-frequency cut-off of the bath which will play an important role in 
what follows. 

In a real experiment of the type discussed the preparation pulse takes 
a finite time, t p . For instance, the (-7r/2) x -pulse which transforms the state 
||) — > (ID + ID)> can be accomplished by putting B z = and B x = huj p 
for a time span t p = it/2uj p . During this time the bath oscillators partially 
adjust to the changing spin state. The oscillators with (high) frequencies, 
ujj 3> lo p , follow the spin adiabatically. In contrast, the oscillators with low 
frequency, loj <C w p , do not change their state. Assuming that the oscillators 
can be split into these two groups, we see that just after the (7r/2) a; -pulse the 

state of the system is "^75 (iD ® g|) + ID ® Sj 1 )) ® gj-) where the super- 
scripts 'h' and T refer to high- and low-frequency oscillators, respectively. 
Thus, we arrive at an initial state with only the low-frequency oscillators 
factorized from the spin. For the off-diagonal element of the density matrix 
we obtain 

\(* + {t))\ = Z(uc,u P )\P Up (t)\, (17) 
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where Z(uj c ,u> p ) = \(g^\g^)\ and Pu> p (t) is given by the same expressions as 
before, except that the high-frequency cut-off is reduced to uj p . 

The high frequency oscillators still contribute to the reduction of | 
- via the factor Z(lu c , uj p ) - however, this effect is reversible. To illustrate 
this we consider the following continuation of the experiment. After the 
(7r/2) pulse we allow for a free evolution of the system during time t with 
magnetic field B z = AE along z axis. Then we apply a (— vr/2) pulse and 
measure a z . Without dissipation the result would be (a z ) = cos(AEt). 
With dissipation the state of the system after time t is 



iAEt/2 



V2 



IT) ® 



|i>®|gf)®e- <w i*/ R |g 1 T )) . 

(18) 

After the (— ir/2) pulse (also of width it/2uj p ) we obtain the following state: 



iAEt/2 



^It>®|gf)®( 



g l\ +e -iAEt/2 e -iH l t/h 



+ 



gl\ + e -iABt/2 e -iW i t/?l 



(19) 



From this we finally get (a z ) = Re P Wp (i)e 



iAEt 



Thus the amplitude 



of the coherent oscillations of a z is reduced only by the factor \Puj p (t)\ 
associated with slow oscillators. The high frequency factor Z(uj c ,uj p ) does 
not appear. To interpret this result we note that we could have discussed the 
experiment in terms of renormalized spins |f) = | Dlgj 1 ) and ||) = | Dl^j 1 ), 
and assuming that the high frequency cutoff of the bath is u) p . 

It is interesting to compare further the two scenarios with instantaneous 
and finite-time preparation further. The time evolution after the instanta- 
neous preparation is governed by P Uc (t). For T <C uj p and t S> l/co p and 
arbitrary spectral density J{ui) it satisfies the following relation: Pu> c (t) = 



Z 2 (uj c ,u p )P u)p (t), which follows from exp(-iHit/h) g\) |<yf|gf)| 2 
for t S> l/ojp. Thus for the instantaneous preparation the reduction due 
to the high frequency oscillators is equal to Z 2 (lo c ,lo p ), while a look at the 
finite-time preparation result (17) shows that in this case the reduction is 
weaker, given by a single power of Z(lu c ,uj p ) only. Moreover, in the slow 
preparation experiment the factor Z(tu c , 

two 'simple' wave functions, gi 1 ) and 



p ) originates from the overlap of 
gj 1 V which can be further adiabati- 



cally manipulated, as described above, and this reduction can be recovered. 
This effect is to be interpreted as a renormalization. On the other hand, 
for the instantaneous preparation the high frequency contribution to the 
dephasing originates from the overlap of the states giMand e -iW J-*/ R 



The latter is a complicated excited state of the bath with many nonzero 
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amplitudes evolving with different frequencies. There is no simple (macro- 
scopic) way to reverse the dephasing associated with this state. Thus we 
observe that the time scale of the manipulating pulses determines the bor- 
der between the oscillators responsible for dephasing and the oscillators 
responsible for renormalization. 



6. Response functions 

In the limit 6 = we can also calculate exactly the linear response of 
t~x = &x to a weak magnetic field in the x-direction, B x (t): 



Xxx(t) = j: Q(t)(T x (t)T x (0) - T x (0)T x (t)) . (20) 

Using the equilibrium density matrix 

P eq = (1 + e-^r 1 [|T) (T| ® p T + e-f> AE |I) (I| ® ^] , (21) 

where p-^ oc exp(— /3Ti.-\) is the bath density matrix adjusted to the spin 
state |t), and similar for pi, we obtain the susceptibility 



Xxx(t) 



2?r 1 e(t) 

1 + e -fSAE 



Im 



(22) 



Thus, to calculate the response function one has to use the full factor P LOc (t), 
corresponding to the instantaneous preparation of an initial state. This can 
be understood by looking at the Kubo formula (20). The operator r x (0) 
flips only the bare spin without touching the oscillators, as if an infinitely 
sharp (n/2) pulse was applied. 

The imaginary part of the Fourier transform of x(*)> which describes 
dissipation, is 

*"» = 2(1 + 1-pae) [ P (^ " AE ^ + ^ AE P(^ + AE)] - . 

(23) 

At T = and positive values of u we use the expression for P(E) from 
Ref. [1] to obtain 



2a- 1 



XLH = ^P(Jiu - AE) = @(hco - AE) 2 ^ (hu - AE) 

(24) 

We observe that the dissipative part x'Lx has a gap AE, which corresponds 
to the minimal energy needed to flip the spin, and a power-law behavior 
as uj approaches the threshold. This behavior of x'xxi 00 ) * s known from the 
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orthogonality catastrophe scenario [5]. It implies that the ground state of 
the oscillator bath for different spin states, |g-f) and |gj), are macroscopically 
orthogonal. In particular, for an Ohmic bath we recover the behavior typical 
for the problem of X-ray absorption in metals [5] . 

As characterizes the dissipation in the system (absorption of en- 

ergy from the perturbing magnetic field) it is interesting to understand the 
respective roles of high- and low-frequency oscillators. We use the spectral 
decomposition for x" at T = 0, 

X » = ttJ2 I <0| r x \n) | 2 [S(u - E n ) - 5{u + E n )} , (25) 

n 

where n denotes exact eigenstates of the system. These are ||) |n-f) and 
||) \ni), where |nj) and |nj) denote the excited (multi-oscillator) states of 
the Hamiltonians H\ and 7i[. The ground state is ||) \g-\) and the only 
excited states that contribute to x"(w) are ||) with Tt^ |n^) = [u — 
AE) \ni). This means that all the oscillators with frequencies Uj > u — AE 
have to be in the ground state. Therefore we obtain for u> p > oj — AE 

x^H = ^ 2 (^,^)^ p H • (26) 

To interpret this result we generalize the coupling to the magnetic field 
by introducing a g factor: H ult = —(g/2)5B x (t)o~ x . Then, if the applied 
magnetic field can be independently measured, the observable quantity cor- 
responding, e.g., to the energy absorption is 

Xxx(t) = g 2 ^e(t)([a x (t),a x (0)}) . (27) 

Thus, Eq. (26) tells us that by measuring the response of the spin at fre- 
quencies to < lo p + AE we cannot distinguish between a model with upper 
cutoff uj c and g = 1 and a model with cutoff lo p and g = Z(lo c , uj p ). This 
is the usual situation in the renormalization group context. Thus, again, 
we note that high-frequency oscillators are naturally associated with renor- 
malization effects. 

We collect the results for x" f° r various spectra in Table 2. The results 
are shown for temperatures lower and higher than the crossover tempera- 
ture introduced in Section 4. 



7. Summary 

The examples presented above show that for a quantum two-state system 
with a non-equilibrium initial state, described by a factorized initial den- 
sity matrix, dephasing persists down to zero bath temperature. An Ohmic 
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TABLE 2. Response functions for different bath spectra. 



Bath spectrum 

Tf \ 7T 1 — S S 

J {iO) — 2QLU UJ 

x9(u c — u) 



Ohmic 
s = f 



Sub-Ohmic 
< s < 1 



Super-Ohmic 
f < s < 2 



Response function 

W*) = *©(t)<M*Wo)]> 



x"x(u), T > T* 




« A) 

ctT 

(aT) 2 +(u-AE) 2 




AE 

ct^p 8(^-Afi) 

(aL; )2 + ( 1J -A£;- Q ^ ln(...)) 2 eX P 



(u-A£-au ln(...)) z 




environment leads to a power-law dephasing at T = 0, while a sub-Ohmic 
bath yields exponential dephasing. The reason is that the factorized initial 
state, even with the bath in the ground state of the bath Hamiltonian, is 
actually a superposition of many excited states of the total coupled system. 
In a real experiment only a part of the environment, the oscillators with low 
frequencies, can be prepared factorized from the two level system. These os- 
cillators still lead to dephasing, whereas the high-frequency oscillators lead 
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to renormalization effects. The examples demonstrate that experimental 
conditions, e.g., details of the system's initial state preparation, determine 
which part of the environment contributes to dephasing and which part 
leads to renormalization. The finite preparation time ~ l/co p also intro- 
duces a natural high-frequency cutoff in the description of dephasing effects. 
We have further demonstrated that dephasing and renormalization effects 
influence the response functions of the two level system. We noted that 
they exhibit features known for the orthogonality catastrophe, including a 
power-law divergence above a threshold. 
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